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ABSTRACT 

PSR B1259— 63 is a middle-aged radio pulsar {P ~ 48 ms, r 



330 kyr, £; = 8.3 x lO^^ erg 
0.87) with a high-mass Be companion, SS 2883. 



s ) in an eccentric binary (Porb = 3.4 yr, e 



We observed the binary near apastron with the Chandra ACIS detector on 2009 May 14 for 28 
ks. In addition to the previously studied pointlike source at the pulsar's position, we detected 
extended emission on the south-southwest side of this source. The pointlike source spectrum can 
be described by the absorbed power-law model with the hydrogen column density A^h = (2.5 ± 
0.6) X 10^-^cm~^, photon index F — 1.6 ±0.1, and luminosity Lo.5-8kcV ~ 1-3 x 10^"^ ergs s~^, 
where d^ is the distance scaled to 3 kpc. This emission likely includes an unresolved part of the 
pulsar wind nebula (PWN) created by the colliding winds from the pulsar and the Be companion, 
and a contribution from the pulsar magnetosphere. The extended emission apparently consists 
of two components. The highly significant compact component looks like a southward extension 
of the pointlike source image, seen up to ~ 4" from the pulsar position. Its spectrum has about 
the same slope as the pointlike source spectrum, while its luminosity is a factor of 10 lower. 
We also detected an elongated feature extended ^ 15" southwest of the pulsar, but significance 
of this detection is marginal. We tentatively interpret the resolved compact PWN component 
as a shocked pulsar wind blown out of the binary by the wind of the Be component, while the 
elongated component could be a pulsar jet. 

Subject headings: pulsars: individual (PSR B1259-63) - stars: neutron - X-rays: binaries 



1. Introduction 

The outstanding resolution of Chandra has 
made possible detailed investigations of pulsar 
wind nebulae (PWNe). PWNe are formed by the 
relativistic pulsar wind (PW), which shocks in 
the ambient medium and emits synchrotron ra- 
diation. The PWN morphology depends on the 
Mach number, A4 = vjcs-, where v is the pulsar's 
velocity relative to the ambient medium, and Cs 
is the speed of sound. PWNe created by subson- 
ically moving pulsars, such as the famous Crab, 



exhibit "torus-jet" morphology, where the torus 
represents a shocked equatorial wind beyond the 
ring- like termination shock (TS), while the jets 
along the pulsar's spin axis are polar outflows. 
At ^ 1, when the ram pressure exceeds the 
ambient pressure, Pram = Pamht^^ > iPa.htn, the 
TS acquires a bullet-like sh ape (jGaensler et al 



2004 iBucciantini et al.ll2005[ ), and the PWN ex 
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hibits a cometary bow-shock morphology, with a 
long tail behind the moving pulsar. Spectacu- 
lar examples of such PWNe are created b y the 
pulsars J174 7-2958 (ICaensler et al.l |2004| ) and 
J1509-5850 (|Kargaltsev et al.ll2008l) . The tails 
may coexist with pulsar jets, often aligned with 



the p roper motion direction (e.g., iPavlov et al 



20031 ). 



An interesting modification of bow-shock PWNe 
is expected when the pulsar is in a high-mass bi- 
nary system. In such a system, the pulsar moves 
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in the wind emanating from the high-mass com- 
panion, and the bow-shock PWN is produced by 
the collision of the relativistic PW with the dense 
nonrelativistic wind of the companion (instead of 
the interstellar medium (ISM) in the case of a soli- 
tary pulsar). The overall appearance, luminosity, 
and spectrum of such a PWN depend on the ratio 
of the momentum fluxes carried by the winds, the 
velocities Vw and Vp of the companion's wind and 
the pulsar's orbital motion, the temperature and 
density of the companion's wind, the anisotropy 
of the outflows, and the s eparation between the 
pulsar and the companion (jTavani fc Aronslll997 . 
hereafter TA97). As these parameters vary along 
the orbit, the PWN's luminosity, spectrum, shape, 
size, and orientation should also vary. If such a 
PWN can be resolved, observations of the variable 
PWN morphology would allow a more definitive 
study of the wind properties than observations of 
an unresolved PWN -I- binary system. 

A promising target for such studies is the well- 
known binary system PSR B1259-63/SS2883. 
PSR B1259-63 (hereafter B1259) is a middle- 
aged radio pulsar, with the period P = 47.8 ms, 
spin-down age t = P/2P — 330 kyr, spin-down 
power E = 8.3 x ip ^s g^-^g g-i^ ^jj^j m agnetic field 
B = 3.3 X 10" G (|Wang et al.ll2004l ). It is in an 
eccentric, wide binary (Pbin = 1237 d w 3.4 yr, 
e = 0.87, a sini = 3.9 x lO^^ cm, i « 36°; 
Johnston et al.lll992l |l994l) with a lOth-magnitude 
Be star SS2883 (M^ « lOAf©, i?* « 6Rq; 



Johnston et al.l 119941 ). The distance estimates 
based on the dispersion measure (DM = 147 
pc cm~'^) give d = 4.6 kpc for the model of 
the Galactic el e ctron density distribution by 
Tavlor fc CordesI (|l993l) . and d = 2.75t^ - ^j, kp c 
for the N E2001 model bylCordes fc Lazioi (l2002l). 



However, iJohnston et al. I (1994') place the system 
at cZ w 1.5 kpc, in the Saggitarius arm, based 
on the properties of the optical companion. As 
the distance remains somewhat uncertain, we will 
scale it to 3 kpc below, d — Sd^ kpc. 



Discovered with ROSAT (|Cominskv et al.lll994[ ) 
X-ray e mission from B125 9 has been studied wit h 



ASCA (iKaspj et al.1ll995t Iffiravama et al.l ll999l). 



Beppo-SAX (iNicastro et aL^ Tl998^. INTEGRAL 



Shaw et al."2004'), X MM-Newton jC hernya kova et al 



2006,. ,2009 1. Suzaku (lUchivama et al.. .20091 an d 



Swift 
TeV 



and Chandra (jChernvakova et al. 

from this system near perias- 



20091) . 



tron was detected wi t h the HESS observatory 
()Aharonian et al.ll2005i l2009l ). The multi-year ob- 
servations have shown that its X-ray flux and spec- 
trum vary with orbital phase. For instance, the 
1 — 10 keV energy flux is at minimum, 1.0 x 10~^^ 



erg cm ^ s ^ 



(i.e., Lx = l-lx 10^3^2 



ergs s 



at apastron, while it is a factor of 20 higher nine 
days before periastron. The slope of its power-law 
(PL) spectrum also varies, from F w 1.7 (at apas- 
tron) to 1.2 (25 days before periastron) to 2.0 (at 
periastron). No pulsations with the radio pulsar 
period have been detected from the system out- 
side the radio range. The orbital variations of the 
X-ray flux and spectrum have been explained by a 
model in which the X-ray emission is generated in 
a PWN formed by the colliding winds; the PWN 
becomes more compact and luminous when the 
pulsar crosses the equatorial disk outflow (angu- 
lar half-width of about 20°) of the Be companion, 
tilted with res pect to the orbital plane b y an angle 
~ 25° fTA97: IChernvakova et aLll2006l ). There is 
no general consensus on the X-ray emission mech- 
anism (synchrotron emission or inverse Compton 
(IC) scattering of companion's optical/UV pho- 
tons). 

It would be interesting to study the unusual 
PWN of B1259 in X-rays, but, at realistic wind 
parameters, the characteristic size of the TS sur- 
face (< O'.'Ol) in the B1259 PWN is too small to 
resolve it even with Chandra. On the other hand, 
we know from the models of bow-shock PWNe and 
observations of supersonic solitary pulsars that the 
flow speed in the coUimatcd tail behind the moving 
pulsar can be so high that the length of the observ- 
able part of the tail can exceed the charac teristic 



TS s ize by a very large factor (Kargaltsc v et al 



20081 ) ■ despite the fast synchrotron and IC cool- 



emission 



ing of the emitting relativistic particles. The tail 
in the immediate vicinity of the binary should be 
oriented along the vector v^ — Vp, i.e., almost along 
the line connecting the Be star and the pulsar (as 
Vw ^ Vp). The tail should rotate with the binary 
period, so that its large-scale appearance would 
resemble an Archimedes spiral ( Dubus 20061) . Op- 
timal binary phases to search for the outflow are 
near the apastron, when the pointlike source emis- 
sion is the faintest, and the size of the bow-shock 
PWN is the largest. We should note that the four 
Chan dra observations taken ne ar the periastron of 
2007 (|Chernvakova et al.ll2009l) were not suitable 
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for detecting extended emission because they were 
short (3-7 ks) and used either Continuous Clock- 
ing mode (which produces only one-dimensional 
image) or High Energy Transmission Grating. In 
addition to the tail outflow, a high-resolution im- 
age may show pulsar jets, which can penetrate 
through the companion's wind, and whose orien- 
tation would not show any orbital dependence. 

To search for extended X-ray emission from the 
PWN tail and pulsar jets, and to better under- 
stand the nature of stellar and pulsar winds and 
their interaction, we observed B1259 with Chan- 
dra near apastron of 2009. We describe this ob- 
servation and the data analysis in Section 2, and 
discuss our results in Section 3. 

2. Observations and Data Analysis 

We observed B1259 with the Advanced CCD 
Imaging Spectrometer (ACIS) on board Chandra 
on 2009 May 14 (ObsID 10089), 39 days after apas- 
tron (i.e., at the binary phase — 0.532, counted 
from periastron, or the true anomaly 9 = 181.63°). 
The target was observed for 28.28 ks in timed ex- 
posure mode using "very faint" telemetry format. 
To minimize the pile-up effects, we used a custom 
124 pixel subarray and turned off all the ACIS 
chips but 13. The frame time in this configuration 
is 0.44104 s (0.4 s exposure time and 0.04104 s the 
image transfer time). The Y-offset of — 0.'2 and 
Z-offset of — 0f25 put the target 29" off-axis, at a 
distance of > 20" from the nearest boundary of the 
w 8' X 1' field of view (FOV), with allowance for 
the 8" amplitude dither. There were no significant 
flares during the observation. The useful scientific 
exposure time (livetime) was 25.65 ks. The data 
were reduced and analyzed with the Chandra In- 
teractive Analysis Observations (CIAO) package 
(ver. 4.2), using CALDB 4.2.0. We used Chan- 
dra Ray Tracer (ChaRT^I and MARJ^ software 
for image analysis, and XSPEC (ver. 12.5.1) for 
spectral analysis. 

2.1. Image Analysis 

To produce images of the pulsar and its vicin- 
ity at subpixel resolution, we removed the pipeline 
pixel randomization and applied the subpixel res- 



^ See [http : //cxc . harvard ■ edu/chart/threads/index ■ html [ 
^See |http : //space .mit . edu/CXC/HARX/'] 
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Fig. 1.— ACIS-I3 images of B1259. Top: Original 
(unbinned, unsmoothed) image in the 0.5-8 keV 
band. The 3"-wide strips northwest and south- 
east of the pulsar are used for the analysis of the 
trailed image of the pointlike source. The 36" x 5" 
rectangle northeast and southwest of the pulsar is 
used for the analysis of extended emission. Mid- 
dle: The 0.5-8 keV image binned to a pixel size 
of 0"98 and smoothed with csmooth to emphasize 
the extended emission. Bottom: Unbinned, un- 
smoothed image in the 1-7 keV band; the semicir- 
cles are used for the imaging and spectral analyses 
of extended emission (see text and Table 1). 
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olution tool to s plit-pixel events in the image 
(|Mori et al.ll200ll : iTsunemi et al.ll200ll ). Figure [D 



shows images of the region around B1259 on the 13 
chip in the 0.5-8 keV and 1-7 keV energy bands. 
The coordinates of the bright compact source at 
the image center (2050 counts in the 1'.'5 radius 
aperture in the 0.5-8 keV band) are R.A.(J2000) 
= 13h02™4??625, Decl.(J2000) = -63°50'08'.'38, as 
measured by the CIAO wavdetect tool. This po- 
sition differs by 0'.'36 from the radio pulsar posi- 
tion, R.A.(J2000) = 13'^02'"47?65(1), Decl.(J2000) 
= -63°50^08y7(l) , for the epoch MJD 50,357 
( Wang et al.ll20o3 ): such a difference is within the 
uncertainty of the absolute Chandra astrometry 
(e.g., 0'.'62 at the 90% confidence level for sources 
within 2' from the optical axis on the ACIS-I3 
chili). 

In the image we see two extended structures 
apparently associated with the pulsar: a narrow 
linear structure seen on both sides of the pul- 
sar along the northwest-southeast direction, and 
a one-sided south-southwest extension of point- 
like source image. As the linear structure is ori- 
ented along the CCD column direction, it could 
be naturally interpreted as the "trailed image" of 
the pointlike source, formed by the source pho- 
tons arriving during the CCD charge transfer to 
the framestore region. To confirm this explana- 
tion and check that no additional extended emis- 
sion is required to produce the observed structure, 
we measured the number of counts in a 3"-wide 
strip, which goes parallel to the chip columns from 
one edge of the subarray to the other (excluding 
a 7" gap with the pulsar in the middle; see the 
upper panel of Figure [T|), and compared it with 
the number of counts expected in the trailed im- 
age. The background surface brightness, 0.044 
counts arcsec"^ in the 0.5-8 keV band, was es- 
timated in a 25" radius circle located about 50" 
northeast from the pulsar. We found 44 counts 
in the strip, of which 9.2 counts belong to the 
background, i.e., there are 34.8 ± 6.7 background- 
subtracted counts. The number of counts ex- 
pected in the trailed image can be estimated as 

A^trail = C'srcitrail(?T-strip/nchip) = 22.5 ± 0.6, whcrC 

C'src = 79.9 ±1.8 counts ks^^ is the source count 
rate in the 1"5 aperture, ttraii — 2.63 ks is the 



^See Figure 5.4 in the Chandra Pro- 
poser's Observatory Guide (POG), ver. 12 



|http : 1 1 cue . harvard ■ edu/proposer/POG/ [ 



total image transfer time during the observation, 
JT-strip = 110 pixels is the strip length in chip co- 
ordinates, and richip = 1024 pixels is the chip size. 
Although the number of counts detected in the 
strip exceeds A'^traii by 12.3±6.7 counts, the statis- 
tical significance of the excess is low (1.8cr), which 
means that the northwest-southeast linear struc- 
ture is indeed the trailed image, without a signifi- 
cant contribution from an extended source. 

To evaluate the statistical significance of the 
other extended feature (south-southwest of the 
pulsar, see Figure 1), we simulated an observa- 
tion of a point source with the same properties as 
the observed ones using the ChaRT and MARX 
packages. To reduce the statistical errors of the 
simulation, we increased the exposure time to 200 
ks when running ChaRT and then rescaled the 
image to compare with the data. We have also 
tried several values of the DitherBlur parameter 
of MARX to match the simulated point spread 
function (PSF) with the real data and found the 
best match at the value of 0'.'2. 

The extended emission apparently consists of a 
compact part, which is seen up to about 4" south 
of the pulsar, and a narrow elongated structure 
of about 10"-15" length southwest of the pulsar, 
which might bend northward at the apparent end. 
To estimate the significance of these structures, we 
investigated the count distributions in two semi- 
circles of 5" radius (north and south of the pul- 
sar; see the bottom panel of Figure [T]) and in 
the 36" X 5" rectangular region (oriented in the 
northeast-southwest direction, i.e., along the elon- 
gated structure; see the top panel of Figure [ij. 

Figure [5] shows the observed and simulated ra- 
dial distributions of counts in the northern and 
southern semicircles, in the 1-7 keV band. We 
see that in the northern semicircle the observed 
distribution is rather close to the simulated one, 
while in the southern semicircle there is a signifi- 
cant excess of the observed counts with respect to 
both the simulated distribution and the observed 
counts in the northern semicircle. For instance, for 
radii between 0'.'75 and 4", the numbers of counts 
in the southern and northern semicircles are 276 
and 160, respectively, i.e., the excess is 116 ± 21 
counts, significant at the 5.6cr level. 

Figure [3] shows the distribution of 0.5-8 keV 
counts along the length I of the 36" x 5" rect- 
angle (I — ^ corresponds to the pulsar position; 
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Fig. 2. — Observed radial count distributions 
(histograms) in the southern (top) and northern 
(bottom) semicircles (see the bottom panel in Fig- 
ure [ij, in the 1-7 keV band, are compared with 
the ChaRT + MARX simulation of a point source 
(dashed line). 



I > southwest of the pulsar). We see an excess 
of observed counts over the simulated profile in 
the southwest half of the rectangle. For instance, 
at 1" < I < 18" there are 158 counts; subtract- 
ing the contribution from the PSF tail and back- 
ground (65.2 and 3.7 counts, respectively), we ob- 
tain the excess of 89 ± 13 counts, significant at 
the 6.9(7 level (if only statistical uncertainties are 
taken into account). The difference between the 
numbers of observed counts in the rectangle seg- 
ments symmetric with respect to the pulsar posi- 
tion confirms the excess. For instance, the differ- 
ence of 158 — 91 = 67 ± 16 counts between the 
1" <l < 18" and -18" < I < -1" is significant at 
the 4(7 level. 

The count excess southwest of the pulsar is. 
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Fig. 3. — Observed count distribution (histogram) 
along the length of the 36" x 5" rectangle shown 
in the top panel of Figure [T] (approximately in the 
northeast-southwest direction) is compared with 
the ChaRT -f MARX simulation of a point source 
(dashed line). The positive and negative distances 
correspond to the southwest and northeast direc- 
tions. The horizontal dashed line shows the back- 
ground level (0.044 counts arcsec^^), while the 
dash-dot line corresponds to the sum of the back- 
ground and the simulation. 



however, mostly due to the brightening in the 
immediate pulsar vicinity, at I w l"-2", while 
the statistical significance of the apparent excess 
at I > 5" (the alleged elongated component) is 
marginal. For instance, in the 1-7 keV band, we 
detected 15 and 3 counts in the 5" < I < 18" and 
— 18" < I < —5" segments, respectively, so that 
the difference, 12.0 ± 4.2 counts, is significant at 
the 2.9(7 level, too low to firmly rule out a statis- 
tical fluctuation. 

The presence of asymmetric extended emis- 
sion at distances > 1" from the pulsar is also 
supported by Figure 4, which shows the result 
of subtraction of the simulated point source im- 
age from the observed one. In addition to the 
asymmetric extended emission, the subtracted im- 
age shows some azimuthally asymmetric system- 
atic residuals of both signs within 1" from the 
point sour ce position, resembling t hose reported 
recent ly by Juda fc Karovska ( 2010l ) and Kashvap 
(|2010tn . This means that the actual PSF core is 



* We should note that the spatial structure of the residuals 
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different from the azimuthally symmetric model; 
this difference, however, makes no effect beyond 
1" from the point source positiord. 




Fig. 4.— Top: ACIS-I3 image of B1259 with 
the simulated point source contribution sub- 
tracted. The ChaRT+MARX simulated image 
(with DitherBlur = 0'.'2) was centered on the 
wavdetect position (see Section 2.1). Both the ob- 
served and simulated images were rebinned to 1 /8 
of the original ACIS pixel size. The residual struc- 
ture within the 1" radius circle is caused by imper- 
fections of the simulated PSF core (see Kashyap 
2010). Bottom: The same image smoothed with 
the C.'IQ Gaussian kernel. 



In principle, an asymmetry of a point source 
image can be caused by the tilt of the inner- 
most Chandra telescope mirrors (the mirror pair 
6) with respect to the other niirrors (see Section 
4.2.3 of POG, and lJerius eraLll2004[ ). It has been 
obs erved in bright sou rces (e.g., the Vela pulsar; 



see iPavlov et al.ll2001l ) as an image elongation be- 
tween the —Z and +Y directions (in spacecraft 
coordinates). In our case, the extended emis- 
sion is seen between the -\-Z and —Y directions. 
Moreover, the effect of the innermost mirror mis- 
alignment is significant for higher photon ener- 
gies {E > 5 keV), while the asymmetry we ob- 
serve in the B1259 image is about the same in 
the 0.5-1.75 and 1.75-8.0 keV images (which have 
approximately the same numbers of counts), and 
our spectrum has too few high-energy counts that 
could cause this asymmetry. Therefore, the ob- 
served one-sided extended emission is not an arte- 
fact caused by the mirror misalignment. We sug- 
gest that this emission can be interpreted as an 
extended PWN of B1259 and conclude that the 
compact component of this PWN, up to ~ 4" from 
the pulsar, is firmly detected, while a longer ob- 
servation is needed to confirm the reality of the 
elongated component. 
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Fig. 5. — Spectrum of the pointlike source fitted 
with the absorbed PL model in the energy range 
of 0.5-8 keV. 



in Figure 4 differs significantly from that reported by these 
authors, but the structure strongly depends on the choice 
of the image center within the ACIS pixel. 
^See |http: //cxc. harvard ■ edu/ciao/caveats/psf .artifact .html#figT| 
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2.2. Spectral Analysis 

2.2.1. Pointlike source 

For the spectral analysis of the pointlike source 
(which likely includes the unresolved part of the 
PWN, in addition to the pulsar), we extracted 
1874 counts in the 0.5-8 keV band from a small cir- 
cular aperture of 1" radius. As the encircled count 
fraction is 88% for such an aperture and energy 
band (according to our ChaRT/MARX simula- 
tions), the aperture-corrected count rate is 83.9 ± 
1.9 counts ks^^, or 0.034 counts per frame, which 
corresponds to a negligibly small pile-up fraction 
(about 1%, according to Chandra PIMM^. The 
background contribution (0.14 counts) is negligi- 
ble in such a small aperture. We binned the spec- 
trum with minimum of 15 counts per bin and fit it 
with the absorbed PL model (see Figure 5). The 
fit is good (x^ = 0.99 for 99 degrees of freedom); 
its results are reported in Table [I] while Figure 6 
shows the confidence contours for the fitting pa- 
rameters. The fit results are consistent with those 
obtained in an XMM-Newton observation near the 
2002 apastron (iVH,2i = 2.9 ± 0.2, F = 1.69 ± 0.04, 
F = (1.04 ± 0.49) x lQ-^^ erg cm"^ s~^ in the 
1-10 keV bancJ3; see IChern vakova et al. ''2006') as 
well as with the ROSAT (C ominsky et al. 1994]) 
and ASCA (jffiravama et al.lll999( ) near-apastron 
observations in 1992 and 1995, respectively. 

2.2.2. Extended emission 

It is difficult to measure the spectrum of the 
alleged PWN not only because of the small num- 
ber of counts but also because the PWN emission 
region overlaps with the wings of the pulsar PSF. 
We, however, attempted to estimate the spectral 
slope and the flux of the PWN emission in the 
two regions shown in the bottom and top panels of 
Figure HJ First, using our ChaRT/MARX simula- 
tions, we extracted the PSF wing spectra in these 
regions and fit each of them with the absorbed PL 
model, freezing the hydrogen column density at 
the best-fit value obtained for the pointlike spec- 
trum, A'h,21 = 2.5 (see Table [1] for the fitting pa- 
rameters). These spectra are slightly harder than 



®See http : //cxc .harvard. edu/toolkit/pimms . jsp 
^Because of the relatively low angular resolution of XMM- 

Newton, those data include both the pointlike source and 

the extended emission. 
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Fig. 6.— Confidence contours (68%, 90%, 99%) 
in the F-iVn {top) and T-N {bottom) planes for the 
ACIS spectrum of the pointlike source, computed 
for the absorbed PL model. The PL normalization 
TV is in units of 10^^ photons cm^^ s~^ keV~^ at 
1 keV. The dashed lines are the lines of constant 
unabsorbed flux in the 0.5-8 keV band (the num- 
bers near the lines are the flux values in units of 



10- 



erg cm 



the PSF core spectrum because the PSF broad- 
ens with increasing energy. Extracting the back- 
ground from the same regions that were used for 
the image analysis, we then fit the spectra in the 
two regions by the sums of two absorbed PL mod- 
els (PSF wings + PWN), with Nu,2i = 2.5 and 
the PSF wings model parameters frozen at the 
previously obtained best-fit values. We see from 
Table [T] that the spectra of the two regions have 
about the same slopes, Fpwn ~ 1-6, statistically 
indistinguishable from the slope of the pointlike 
source, while the PWN flux in the southern semi- 
annulus (S) is a factor of 3 higher than that in 
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Energy (keV) 

Fig. 7. — Spectra of the background-subtracted 
emission in the SW (top) and S (bottom) re- 
gions (see Table 1 and Figure 1). The red his- 
tograms show fits with the absorbed PL models of 
the corresponding PSF wing spectra obtained in 
the ChaRT-l-MARX simulations. The black his- 
tograms show fits of the observed (PSF+PWN) 
spectra with the sum of two absorbed PL models, 
one of which (with fixed parameters) describes the 
contribution of the PSF wings. 

the southwest rectangle (SW). The total flux and 
the luminosity of the alleged PWN are a factor of 
~ 10 lower than those of the pointlike source. 

3. Discussion 

3.1. Pointlike source 

The pointlike source emission can include con- 
tributions from the pulsar, the Be companion, and 
the unresolved part of the PWN created by the col- 
liding winds from the pulsar and the Be compan- 
ion, including the ultracompact intrabinary PWN, 



with a size smaller than the separation between 
the binary components, < Sdg ^ milliarcseconds 
at apastron. Coronal emission from the Be com- 
panion, with an optically thin, thermal plasma 
spectrum corresponding to a temperature of a 
few million kelvins, might give some contribution 
in low-energy channels, but, given the goodness 
of the single-component PL fit, the contribution 
of the coronal emission can be neglected. The 
pulsar emission could include a thermal compo- 
nent, emitted from the neutron star surface, and a 
nonthermal component emitted from the pulsar's 
magnetosphere. Since the observed spectrum ex- 
cellently fits the single PL model, we conclude that 
the pulsar's thermal emission is too soft and faint 
to be seen at the relatively large A'h ~ 2.5 x 10^^ 
cm~^. For instance, the thermal emission from the 
Geminga pulsar, whose spin-down age is similar to 
that of B 1259 and whose bulk surface temperatur e 



< 5 X 10^ K (see, e.g.. iKargaltsev et al.|[2005l) . 



is 

would be virtually undetectable at energies above 
0.5 keV if Geminga were at the distance of B1259. 

For a typical X-ray efficiency of the magneto- 
sphe ric emission, ryx.psr = Lx.os r/E ~ 10~^ (see, 
e.g., IKargaltsev fc Pavlov l2008h . one would ex- 
pect Lx,psr ~ 10'^^ erg s~-^, which is an order of 
magnitude lower than the detected Lo.5-8keV = 
1.3 X 10'^^d§ erg s^^. However, as ?7x,psr shows 
a large scatter (from < 10~^ to ^ 10~^ for 
pulsars with similar valu es of E and t — see 
Kargaltsev fc Pavlovll2009l) . the relatively high X- 



ray luminosity of B1259 does not rule out a sub- 
stantial contribution from the pulsar, especially 
taking into account that the distance may be 
smaller than the assumed 3 kpc. As the magneto- 
spheric emission is expected to be strongly pulsed, 
its contribution could be estimated by measuring 
X-ray pulsati ons near apastron . Such an attempt 
was made by iHiravama et al. ( 19991 ) in observa- 
tions with the ASCA GIS detector, but the up- 
per limit on the pulsed fraction, Pf 40%, was 
not constraining. To conclude, the pointlike emis- 
sion likely includes contributions from the unre- 
solved part of the PWN and the pulsar's magneto- 
sphere, with the total X-ray efficiency fyo.s-skcV = 
1.6 X 10"3(i§. Timing observation with XMM- 
Newton or Chandra could help separate their con- 
tributions. 
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3.2. Extended emission 

Thanks to the unprecedented angular resolu- 
tion of Chandra, our observation of B1259 has re- 
vealed, for the first time, extended X-ray emission 
associated with this remarkable binary. As the X- 
ray emission with a PL spectrum is usually gener- 
ated by relativistic particles, and the pulsar is the 
only natural source of such particles in B1259, we 
interpret the extended emission as a PWN. 

The detected extended PWN emission appears 
to consist of two components. Firmly detected is 
only the relatively compact component south of 
the pulsar, with a characteristic size of ~ 4" and 
luminosity of ^ lO'^^dg erg s~^. The detection of 
the less luminous elongated component, of ^ 15" 
length, is less certain. 

The elongated component could be interpreted 
as a pulsar jet. Its ~ 15" angular length corre- 
sponds to ~ 7 X 10^^ cm at d = 3 kpc, which is 
comparable with the 5 x 10^*" cm and ^ 2 x 10^^ 
cm projected lengths of the X-ray je ts of the Vela 



and Gemin ga pulsars, respectively (jPavlov et al 



10^^ erg s ^, 



IS 



2003L hOldt ). Its luminosity, 
about 10~^ of the pulsar's spin-down power, ver- 
sus ~ 10^^ for the Geminga and Vela pulsars. To 
check this interpretation, longer Chandra obser- 
vations of B1259 are needed, at different binary 
phases. If the detected elongated feature is indeed 
the pulsar jet, we do not expect substantial varia- 
tions in its luminosity and length, except may be 
for a short interval (~ 40 d) around periastron, 
when accretion of the Be c ompanion wind m ight 
quench the pulsar activity (jWang et al.ll2004l ). 

The compact PWN component could be asso- 
ciated with the shocked PW blown out of the bi- 
nary by the wind of the Be companion. Similar 
to the case of a single pulsar moving in the ISM 
with a supersonic velocity, the collision of the two 
winds should create a bow shock within the bi- 
nary (TA97) and a tail directed out of the bi- 
nary approximately along the line connecting the 
Be companion with the pulsaiH. As the pulsar 
moves along the binary orbit, the tail is expected 
to be bent in the direction opposite to the pulsar 
motion, resembling an Archimedes spiral ( Dubud 



'We assume here that the pulsar's velocity (which varies 
from 12 km s~^ at apastron to 176 km s~^ at periastron) 
is much smaller than the velocity of the Be companion's 
wind near the pulsar (see TA97). 



20061 ). The tail's shape can described by the equa- 
tion 



2tt 



(1 -f ecos6'')" 



(1) 

where R and 9 are the radius and the azimuthal 
angle of a point on the tail's central line in the 
reference frame centered at the Be companion, 
Rp and Op are the radial coordinate and the true 
anomaly of the pulsar, and Vfiow(^') is the flow 
speed that may vary along the tail. The shape of 
the tail near apastron, when |1 -I- cos6'| ^ 1 and 
|1 -|-cos0p| ^ 1, is given by the following equation 

R-Rp[6) « ^^|^^(^.-^) 



= 2.3 X 10 



19 ^4iow 



27r 



■cm, (2) 



where t^iow is the average flow speed, and the 
numerical estimate is given for the B1259 case 
[Op = 1.00907r for our observation). The compact 
PWN in the ACIS image looks like a one-sided 
bulge on (or a "blob" attached to) the pointlike 
source rather than a narrow tail. This is expected 
if the angular size of the visible tail only slightly 
exceeds the size of the PSF. In addition, the tail 
may significantly broaden by interactions with the 
ambient medium. 

The luminosity of the compact component of 
the resolved PWN, Ln.s-sknV 1-2 x lO^^dl erg 



is a factor of 



^0.5~8kcV 

10 lower than the lumi- 



nosity of the pointlike source, which includes the 
PWN and pulsar contributions in an unknown 
proportion. The corresponding X-ray efficiency, 
?7o.5-8kcV = Lo,5^skcv/E = 1.4x 10~'*d§, is at the 
lower end of those observed in the tails of super- 
sonically moving solita ry pulsars (see Table 5 in 
Kargaltsev et al. l2008l ) , but one should take into 
account that a substantial part of the alleged tail 
of B1259 cannot be resolved even with Chandra, 
so the true tail luminosity and the X-ray efficiency 
may be considerably larger than these estimates. 

The observed angular size of the compact 
PWN, ~ 4", which corresponds to the projected 
length of '--^ 2 X lO^^da cm, and the fact that the 
PWN is seen only on one side of the pulsar, put a 
constraint on the average flow speed: Vflow ^ O.lc. 
Such high flow speeds have been obtained in nu- 
merical models for tails of bow-shock PWNe of 
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solitary pulsars (e.g., Vhow '~ 0.1c-0.3c in the nar- 
row central channel of the tail, and ~ 0.8c-0.9c in 
the tail's "sheath", in the models by Bucciantini 
et al. 2005), but these speeds may become lower 
at large distances from the pulsar because of the 
interaction with the ambient medium (this effect 
was not included in the numerical models). Ob- 
servations of tails behind supersonically moving 
solitary pulsars also show rather high flow speeds 
— for instance, Vhow 0.02c for the 6 pc tail 



behind PSR J1509-5850 ( Kargaltsev et al 



20081) . 



Even higher flow speeds, corresponding to bulk 
motion L orentz factors 7hii i k up t o ~ 100, were ob- 
tained bv lBogovalov et al. 1(2008') in hydrodynam- 
ical simulations of flows produced by the collision 
of the relativistic (pulsar) and nonrelativistic (Be 
star) winds in a binary system. The ultrarelativis- 
tic flows are formed when the ratio rj of the mo- 
mentum flux from the pulsar to that from the Be 
star exceeds 1.25 x 10~^, in which case the TS sur- 
face becomes "unclosed" (i.e., the back surface of 
the TS bullet disappears). In the case of B1259, 
r] may vary between 10~^ and ^ 1 along the 
binary orbit, being as large as 0.1-1 around apas- 
tron (when the PW interacts with the so-called 
polar component of the Be star outflow). If the 
confined flow indeed reaches such high speeds, it 
can penetrate farther into the ambient medium, 
which makes the suggested interpretation of the 
obser yed asymmet r ic eni ission even more plausi- 
ble. Dubus et al. ( 2010l) note that the proper- 
ties of the observed X-ray and VHE emission from 
B1259 are not compatible with the strong Doppler 
boosting corresponding to such large 7buik- How- 
ever, it is possible that, while the bulk of the ob- 
served high-energy emission is generated within 
(or very close to) the binary, where the flow is 
only mildly relativistic, the extended emission is 
generated by electrons supplied by the accelerated 
flow that slows down at large distances from the 
binary. 

The interpretation of the compact PWN as a 
shocked PW outside the binary can constrain the 
wind paramenters in various models of B1259. 
The X-ray PWN emission could be due to syn- 
chrotron radiation of ultrarelativistic electrons, 
accelerated at the inner shock front within the bi- 
nary (TA97). To generate X-ray photons, E > 0.5 
keV, the electron energy jmc^ and the magnetic 
field B should be high enough: B'^^ > 10^^ (where 



B is in units of gauss). On the other hand, the 
cooling time, icooi, for the relativistic electrons 
in the flow should be longer than the flow time, 
iflow ~ -R/Vflow ^ 3 x_10'^i?i7Vg~5^ s, where Rn = 
i?/(10"cm), V9.5 = Vflow/(10^-^cms-i). Gener- 
ally, both the synchrotron cooling in the magnetic 
field and IC cooling in the radiation field of the 
Be companion can be important. At the very high 
energies of electrons emitting synchrotron X-rays, 
7 > mc^/e 10^ (where e ~ SkT^s 5 eV is 
the mean photon energy), the IC cooling occurs 
in the Klein-Nishina regime. The IC cooling rate 
in this regime is much lower than in the Thomson 
regime: 7ic,kn/7ic,t ~ {mc^ /"fe)'^\n{'ye/mc^) ~ 
10"'* 77^(ln77 -I- 5), and it is much lower than the 
synchrotron cooling rate beyond the binary orbit. 
Therefore, the cooling time is determined by syn- 
chrotron radiation, tcool ^ isyn ^ 5xlO*7"^i?~^ S, 
and it exceeds the flow time at 7-6^ < 16Vg.5i?]^/. 
Taking into account the condition Bj^ 10^^ de- 
rived above, we obtain 



B < lAxlQ-^vli^Rj^^ G, 
7 > 0.9x 10^K,://'i?}(^ 



(3) 



for the typical magnetic field and Lorentz fac- 
tor. The upper limit on magnetic field is much 
lower than the magnetic field in the ultracompact 
intrabinary PWN (e.g., B - 1 G in the TA97 
model), while the lower limit on Lorentz factor 
is close to 7max of the electron energy spectrum 
in the shocked (uncooled) PW within the binary 
(TA97). To reconcile such high energies and low 
magnetic fields with the intrabinary PWN mod- 
els, one has to assume that the radiating electrons 
are additionally accelerated in the tail, in some 
internal shocks or via magnetic field reconnection 
that transforms the magnetic field energy into the 
kinetic energy of particles. 

As noted by Chernyakova et al. (2006, 2009), 
the TA97 model cannot explain some of the ob- 
served properties of the B1259's radiation (such as 
the hardening of the X-ray spectrum near perias- 



tron). These authors (see also Neronov fc Che rnvakova 
20071) discuss another model (fi rst considered by 
Cher nvakov a &: Illarionovlll999l) . in which the X- 
ray emission is generated by the IC scattering 
of stellar photons off moderately relativistic elec- 
trons with Lorentz factors 7 ^ 10-100 (much 
lower than in the TA97 model). As such scat- 
tering occurs in the Thomson regime, the cool- 
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ing rate is 71c 



-(4/3)((TTc/wc^)7^?7ph, where 



t/ph = L*/(47rci?2) = 1.7 X 10U,,38i?r7 eV 
cm"'^ is the photon energy density at the dis- 
tance R — 10^^i?i7 cm from the Be companion 
(approximately equal to the distance from the bi- 
nary a.t R:$> a), and — lO'^^L^^ag erg s^^ is the 
luminosity of the stellar companion. At realistic 
values of the magnetic field, B < SOOL^^^sRij 
/iG, the photon energy density exceeds the mag- 
netic energy density, so the IC cooling is more 
efficient. Taking into account that 71c oc R~^, 
the electron's energy at distance R is given by the 
equation 



liR) = 70 



1 



To 
7c 



1 



(4) 



where 70 = 7(-Ro), 7c = iirmc'^RoViiow/icrTL^) « 
3.7 X 10'^Ro^uV9.5L~lg. If we assume 70 ~ 10-100 
close to the binary orbit at apastron {Rq ^ lO^'' 
cm), then7o/7c (3~30) x lO^^y-^i^ ^j^j^jj^ 

means 

that the total relative energy loss by an elec- 
tron, (70 - 7oo)/7o = (7o/7c)[l + (7o/7c)]"^ « 
7o/7ci is very small at realistic 14iow ^ 10- 
100 km s~^. Moreover, the upper limit on 
the X-ray efficiency of the resolved PWN be- 
tween distances Ri and R2 from the binary, 

[7(^1) - 7(i?2)]/7o - (7o/7c)i?o(i?r' - ^2"') ~ 
10-'^(7o/30)K,75^ (Rll^ - R2.I7), is much smaller 
than the observed value, ^ 10"'''. This means 
that the observed extended X-ray emission can- 
not be explained by the IC scattering, and the 
synchrotron emission of electrons with TeV ener- 
gies is a more plausible interpretation. 

The suggested interpretation of the compact 
PWN as a bent tail of the shocked PW blown 
out of the binary by the wind of the Be compan- 
ion implies that the compact PWN rotates around 
the pointlike source with the binary period, and 
its size and luminosity should depend on binary 
phase. This can be verified by observing this sys- 
tem with Chandra at different binary phases. 

If these observations show the compact PWN at 
the same position with respect to the pulsar, an- 
other interpretation would be required to explain 
the extended emission. In principle, it might be 
a tail or trail of shocked PW behind the moving 
binary blown out by the head-on wind of the am- 
bient ISM matter. The proper motion of B1259 
has not been measured yet, but the speed of the 
high-mass binary should be much lower than that 



of a solitary pulsar, so that it would be hard to 
explain even so short tail/trail as caused by the 
motion of B1259 in the ISM. 

Finally, we note that strong evidence of ex- 
tended emission around another 7-ray binary, the 
microquasar LS 5039, has be en recently found i n 
a Chandra ACIS observation ( Durant et aLllioil ). 
It shows that the compact object in LS 5039 is 
likely a rotation-powered pulsar rather than an ac- 
creting neutron star or a black hole. It also means 
that B1259 is not the only high-mass X-ray binary 
surrounded by an extended PWN. 
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Table 1: Properties of X-ray Emission from B1259 



Region'* 


A^H.21 


r 






Tjiunabsd 
^-14 


-^32° 


Pointlike 


2 'sn+"-*^" 

^•^'^-0.54 


-'^•'-'^-0.11 


171+^^ 


83.0 ±4.5 


1191!^ 


13 


PSF S 

PWN S 


[2.5] 


1 07+U.I8 

-0.17 
1 r7+0.29 
1.0(_g 27 


10±2 

17±4 


9± 1 


11±2 


1.2 


PSF SW 
PWN SW 


[2.5] 


1 34+'J-it> 

i-J'i-O.lS 

l.DO_Q 49 


1+0.8 
^•-'■-0.7 
7 q+3.4 
'•"5-2.7 


3±1 


4±1 


0.4 



Note. — The errors shown represent 90% confidence intervals. The fluxes and hiiniriosities are for the 0.5—8 keV band. 
"The region marked S is the semiannulus 0'.'75 < r < 4" in the semicircle south of the pointlike source, shown in the bottom 
page of Figure 1. The region marked SW is the southwest part, 1" < I < 18", of the 36" x 5" rectangle shown in the upper 
panel of Figure 1. 

''PL normalization in units of 10~® photons keV~^ cm~^ s~^ at 1 keV. 

■^Observed (absorbed) flux in units of 10"^'' erg cm~^ s~^ (in 1" radius aperture for the pointlike source, which contains 86% 
energy fraction). 

''Unabsorbed flux in units of 10^^^ erg cm~^ s~^ (aperture-corrected for the pointlike source). 
^ Luminosity in units of lO^'^dg erg s~^. 
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